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Abstract: Based on the modern representation of the water softening process by the ion exchange method, 
modeling was carried out, an algorithm was developed and the ionite regeneration process was 
calculated in a water softening plant with a fixed ionite layer. To solve the system of equations of the 
model, special, but simple enough for use in engineering calculations, techniques and methods have been 
developed and applied. The constructed model and its solution algorithm can be used to predict the 
process of ionite regeneration in an ion exchange apparatus. It can also be used to calculate the time of 
the ionite regeneration process. 


Keywords: modeling, sorption, regeneration, water softening, ion exchange method, ionite layer, clarified 
water, cationite, anionite, water hardness, sorbent height, kinetic and equilibrium dependencies, 
overshoot concentration. 


Introduction 


The conditions for stopping the regeneration of the sodium-cationite filter depend on many factors, 
such as the concentration of ions in water, the concentration of ions in the ion exchange layer, the length 
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of the ion exchange layer, and others. Usually, the ion exchange layer can be regenerated when the length 
of the ion exchange layer reaches a certain limit, which depends on the concentration of ions in the water 
(the hardness of the filtrate reaches 50% of the hardness of the source water) [1]. Also, some ion- 
exchange materials have the ability to autoregenerate when properly installed and configured. 

The technology of cationite regeneration consists of the following operations: loosening, 
regeneration and washing. 

Cationite regeneration is carried out to restore the ion-exchange properties of the ion-exchange 
mixture. Cationite is an ion-exchange agent that absorbs and removes impurities from water, such as 
hardness, salts and other undesirable substances [2]. During the ion exchange process, the activity of the 
cationite decreases, which leads to the loss of its ion exchange properties. Therefore, in order to maintain 
the ion-exchange properties of cationite, it is necessary to periodically regenerate it. Ionite is regenerated 
with 8% sodium chloride solutions according to a direct-flow scheme: 

CaR2+2NaCl=2NaR+CaClh 

Cationite regeneration usually occurs when its activity falls below a certain level [3]. This usually 
happens after a sufficient amount of impurities has collected in the water, which the cationite can no 
longer remove. Some systems have set activity limits that, when triggered, require cationite regeneration. 
In other systems, regeneration is performed according to a specific schedule. For example, weekly or 
monthly. 

The washing of the ionite material is carried out in order to clean it from the impurities present on it. 
Washing allows you to remove salts and other harmful impurities from the ionite material that accumulate 
on the surface of the ionite material during the ion exchange process. Washing helps to keep the ionite 
material in good condition and improve its performance. It also helps to prevent clogging of the ionite 
material and improve the quality of the treated water [4]. 

Let's make a mathematical description of the ionite regeneration process in the apparatus, the 
scheme of which is shown in Fig.1. 
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Fig. 1. The scheme of operation of the device at the stage of regeneration: 
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1 - housing; 2 - partition; 3-cover; 4,5,13- fittings; 6,9 - upper and lower switchgear; 7,8,10,11,12 — 
chambers; 


A layer of cationite is located in chamber 11. Cationite is a spherical grain with a diameter of 2ro. A 
layer of anthracite is loaded into chamber 8. During ionite regeneration, exchange tanks are used to 
preserve the ionite's ionic properties and increase system performance. During the regeneration of this 
type of filters, the regenerating solution is fed through the fittings 4 and 5 through chambers 12 and 7 into 
chambers 8 and 11, passes through a layer of strong acid cationite (CU-2-8) and anthracite into chamber 
9, and is removed through the fitting 13. 

After regeneration, the process of flushing through the fittings 4 and 5 follows in the same order as 
during regeneration. 

The model of the regeneration process was based on the following assumptions: 

ehe speed of the ion exchange process is limited by both external and internal diffusion; 

ehe one-dimensional flow of the liquid phase in the apparatus depends only on the x coordinate; 

ethe change in the concentration of the solution is a function of the average velocity of the liquid 
phase, only the longitudinal mixing and mass transfer process between the ionite and the solution, the 
equilibrium of which is described by the equation of the linear Henry isotherm; 

The mathematical description of the process, based on the accepted assumptions, consists of the 
following equations: 


material balance: 
mn 2 
rciean | aPS che yee ge . (1) 
OT Of Ox Ox 


an equation reflecting the kinetics of diffusion into a spherical particle: 
OC = OC 20C 
ee a ee 2D, 
OT oan r ar ) 
the equations of the nonlinear Henry isotherm: 


C,=1C (3) 
an equation reflecting the relationship of the local concentration of a substance in the solid phase 
with its average value 


C,., (4) = = i) r°C(t,r)dr (4) 
'% 0 
The mathematical description includes initial and boundary conditions: 
C(z,x)|__, =Co (5) 
C(t, x)| __, =Cin (6) 
OC(t, x) _0 (7) 
Ox x=H 
C(t.r)|_, = Co) (8) 
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eC(r,r)} _ B 
or r= - Dy 7 a 
6C (tr, r) ~0; (10) 
or r=0 


To solve this mathematical description, as in the case of the sorption process, the ionite layer along 
the entire height of the apparatus is divided into elementary layers each with a height of h = H/n. The 
time of the technological process will look like sequentially connected time intervals [5]. 

The second component of equation (1) determines the amount of runoff depending on the ion 
exchange adsorption process, and it can be replaced by an equation in the form of a finite difference: 


dy Ay 
For this equation, the index j determines the number of the time interval, and the difference between 


OCmedi ee Cmedi (tj +hr)—Cmedi (x;) = Kei; fal 1) 


the time intervals is defined as At=1}+1-T). 

Solution of a system of equations (2) — (4) reflecting the kinetics of diffusion into a spherical 
particle, the nonlinear Henry isotherm and an equation reflecting the relationship of the local 
concentration of matter in the solid phase with its average value, taking into account the initial and 
boundary conditions (equations (8) — (10) by analogy with solving the problem of determining the thermal 
conductivity of spherical bodies the forms, provided that the concentration of the target component is 
unevenly distributed, make it possible to determine the amount of runoff in the following form: 


= _ tw. = 
Go~\medi _ 20~"oj F100 B,, exp (- 2 Be) (12) 
0 


ag ag n=—1 
In equation (11), the value of Bz is calculated by the formula 
2 ae 
£2. (Un—Sin fy COS fty,) (13) 


where [Un are the roots of the equation 


n 


1 
te-— mf 


DefG 
The joint solution of the material balance equation with the equations of the initial and boundary 
conditions (5) — (7), taking into account the equation determining the amount of runoff, allows the 
determination of a ratio that allows calculations to be made to determine the spatiotemporal change in the 
concentration of the regenerating solution when it passes a stationary ionite layer. 
In order to facilitate the solution of equation (11), we introduce symbols in the form of 
dimensionless quantities: 


gee Cc =6 = ag — C — ag—G, 
— m No = mm. °. Nea = 0 =. No _ 0 iP 
Cin Cin ag ag 
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D, tT x r 
Fo =F 3z2=5 O<x<h); § =~ (14) 
After entering the above symbols, equation (11) can be rewritten in the following wording: 
ONmed _ Nmed(Fo)-Nmea(Fo+4Fo) = — 
aFo AFo — Kn (15) 
where 
ie 
ef 


As a result of the transformations carried out , the mathematical model of the ionite regeneration 
process can be written as follows: 


an = a2N an 
N(Fo.z 
e Ms = Nin (17) 
No.2) = No (18) 
Fo-0 
aN(Fo,z) = 
i eae (19) 


Based on the above relations, we will solve the nonlinear, inhomogeneous differential equation of 
the material balance (16) with initial and boundary conditions (17) - (19). 
Based on the accepted designations (14), we transform equation (16) to the following form: 
a-e(2) vp 2°N(r2) va Ewce2)) 
c(2.ucer)) + Oe) — raleences) _ rset) (17) 
Based on the Fourier method, we will look for a solution to equation (17) in the form of the 
following series: 


N(z,7) = 2-14, sin( (18) 

Assuming that the function C(t,x) can be decomposed into a Fourier series with respect to the 
variable x, we substitute the relation (19) into the equation (17). 

N(z,r) ?= Din-1 Nn(®D sin (=) (19) 

As a result, we obtain the following relation: 


* me fal 
(> (muco }sin (= )s ad ha - 


n=1 


TENZ 
h 


: 2 
Dpre (Sy (x2) vt (Ey (c2)) 


= SE en (20) 
Defh Defh 
From this equation follows: 
: Zoe 
< N(x) + = 0 (21) 
The solution of this equation in general form can be represented as the following relation: 
_Defr*n*r 
N,@) =_Cle (22) 
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Based on conditions (17) - (19), it is possible to deduce the initial conditions for N(t ) when t=0. 


Then the solution satisfying the initial N (t) =0 conditions of equation (21) can be presented in the 


following form: 
_Defr*n* (zt) 
N,,() =f, te h? dt (23) 
Substituting solutions (23) in a series (18), one can obtain a solution of equation (16) under initial 
conditions (17) - (19) in the following form: 


_Defr'n* (zt) 
AG) = 75 (i Te Pe ar) sin (=) (24) 


If the initial conditions are not equal to zero, it is necessary to add solutions satisfying the initial and 
boundary conditions of the homogeneous equation to solution (24). 
The result can be represented as the following Grin 's function: 
Na(X) *= f, fy GG, &,t,1) d&dz (25) 
As a result, we get the following function: 
Defrn*(z—t) 


aed the solution of oasis (16) under initial and boundary conditions is as follows: 
Defr7n7(r-t} Defrwe® in (r—t} 
sels ei) 


2h? sin(—— )e h 2(—1)?*" sin e 2h? sin(—— 
N.G) = ‘ CF) = 4 Fy 7 / 4 Sy. 4. 
n= Def? n?(t-t) nm Deft’ n= Defrn?(t—t) 
7 2(-1)?*"n?sin(=) 
n=1 Defr=n? 
(27) 


Now let's consider the solution of a nonlinear inhomogeneous differential equation of diffusion 
kinetics (2) with initial and boundary conditions given in the form of equations (8 - (10). 
Using the accepted conditional values, we transform equation (2) to the form: 


= Ne(r, r= ver ( = Nt(z,r) - (28) 


As in the previous case, using the Fourier method, we transform the solution of this problem to the 
form of the Grin function: 


Defm?n?(r-t) 
3 Woo 28 nh? i (=) : (=) 
G(z,zi,t,t) = Din — > sin| >) sins (29) 
In this case, the solution of the diffusion kinetics equation (2) with initial and boundary conditions 
can be written as follows: 


Defr'n(r-t) Defr-n* (rt) 
2n?sin(=)e ” 2(-1)7*"n2sin()e - 2n? sin(***) 
N,,(z, 1) ?= ale Defrn® — n=1 Defren® -( m=1" Derren? )s 
+52, 2(—1)?*"n?sin(=*) 
| Defrin® 
(30) 
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To solve this model, we will use the finite element method - the PDE Toolbox of the Matlab 
application software package. We denote the packages in the PDE toolbox in the following form: 


——— 


t : . . =~ AA 7 


SS he SS ff =3 


Fig.2. Representation of the model layout in the PDE toolbox. 


To determine the interval of possible solutions of the mathematical model (30), we use the 
Neumann method and introduce the corresponding values: 


Fig. 3. The method of determining solutions in the PDE toolbox 


moe “ee od 


——- 


Fig.4. Choosing a solution method in the PDE toolbox 
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To determine the coordinates of the location of possible solutions in the PDE toolbox, it is 
necessary to enter the appropriate model parameters and determine the following dependencies (Fig.4, 
Fig.5): 


— ——— —————— - 7 


| a 


Fig.6. Coordinates of the location of possible solutions in the PDE toolbox 
To determine the range of possible solutions according to the Neumann method, coefficients 
corresponding to the model parameters are introduced: 
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Fig.7. Entering the corresponding coefficients in the PDE toolbox to obtain solutions 
After that, in order to clarify the solutions using the Neumann method, we will construct 
histograms: 


A = wal - | 
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Fig.8. Possible solutions of the model on the PDE toolbox 


14 


Fig.9. Acceptable model solutions on the PDE toolbox 

This means that the permissible solutions of the model are located in the interval (0 — 40) and, as 
can be seen from Figure 9, the closest values of the permissible solutions are in the interval (25 — 35). 

Table 1 shows the values calculated according to the given model and experimental data on 
cationite regeneration. 
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Exchange Na* - Ca** Exchange Na* - Mg** 
Time T, s ree (Time tT, s mi 

experiment calculation experiment calculation 
1 2 3 4 5 6 
240 0,680 0,678 180 0,698 0,721 
300 0,633 (0,629 240 (0,679 0,701 
360 0,574 0,567 300 0,644 0,679 
420 0,492 0,479 360 0,605 (0,633 
480 0,441 0,435 420 0,578 0,538 
540 0,387 (0,369 480 0,546 0,515 
600 0,293 0,283 540 0,502 0,498 
660 0,271 0,268 600 0,467 0,488 
720 0,198 0,189 660 0,413 0,449 
780 0,143 0,147 720 0,387 0,407 
840 0,127 0,138 780 0,341 0,358 
900 0,106 0,109 840 0,305 0,327 
960 0,098 0,101 900 (0,279 0,298 
Table 1 


Experimental data and cationite regeneration values calculated according to the developed model 


The given experimental and calculated data show that the calculated data satisfactorily correspond 
to the experimental ones. The standard deviation of the experimental values from the calculated values 
does not exceed 16%. 

A comparative picture of the two parallel experiments is shown in Figures 10 and 11. 
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Fig.10. Comparison of experimental and calculated values of regeneration curves Ca?*-Na* 
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Fig.11. Comparison of experimental and calculated values of regeneration curves Mg”*-Na* 

Summing up the research, it can be stated that an adequate mathematical model of the ionite 
regeneration process has been developed in a water softening plant with a fixed ionite layer, as eloquently 
evidenced by the proximity of experimental and model-calculated data, as well as the close coincidence of 
the curves shown in Figures 10 and 11. 
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